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Abstract 

The rounding of first order phase transitions by quenched randomness is stated in a form 
which is apphcable to both classical and quantum systems: The free energy, as well as the ground 
state energy, of a spin system on a d-dimensional lattice is continuously differentiable with respect 
to any parameter in the Hamiltonian to which some randomness has been added when d < 2. This 
implies absence of jumps in the associated order parameter, e.g., the magnetization in case of a 
random magnetic field. A similar result applies in cases of continuous symmetry breaking for 
d < 4. Some questions concerning the behavior of related order parameters in such random 
systems are discussed. 

1 Introduction 

The effect of quenched randomness on the equilibrium and transport properties of macroscopic 
systems is a subject of great theoretical and practical interest which is close to Nihat's heart. He 
and his collaborators 1 1 , 2 , 3 1 made important contributions to the study of the changes brought 
about by such randomness in phase transitions occurring in the pure (non-random) system. These 
effects can be profound in low dimensions. Their understanding evolved in a somewhat interesting 
way. As was argued by Imry, Ma and others 1 1] |4] |5l, even a very small amount of randomness 
can change the nature of the transition in dimensions d < 2. The validity of the reasoning involved 
was contentious as it conflicted with the prediction of a dimensional-reduction principle, which was 
arrived at through perturbative calculations carried to all orders [6|. Rigorous work by Imbrie [7J 
and Bricmont and Kupiainen |8| showed that, in the discrete case, the rounding phenomenon does 
not extend beyond two dimensions, as given by the Imry-Ma argument. Aizenman and Wehr ||9l 
[lOJ provided a rigorous proof of the rounding effect in low dimensions, which although limited to 
classical spin systems included a large variety of types of randomness not previously considered. 
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As we have recently shown fTTl, with some modification this proof can be extended to quantum 
systems. Here we present these results in a unified way, and discuss related issues which appear in 
specific examples of interest. 

To formulate the results in a general, but not the most general way, we consider a lattice system 
of generalized spin variables a_ = {ct^^}, x € Z'^, where cr^ is an operator in a finite dimensional 
Hilbert space or, in the classical case, a function taking values in a finite set. The Hamiltonian is of 
the form 

He = IC{g:;J) - h'^K^ - e^rj^n^ = Ti-o - e^rj^Kx (1.1) 

XX X 

where J denotes a set of interaction parameters, {k^} are translates of some local operator kq, and 
h and e are real parameters. The quenched disorder is represented by {rjx}, a family of independent, 
identically distributed random variables with an absolutely continuous distribution with mean zero 
and finite third moment (weaker conditions are possible, but more awkward to state). /C((t; J) may 
be translation invariant and nonrandom, or it can include additional random terms but we will not 
discuss the latter case. It is assumed here that the interactions decay faster than (distance)^'', with 
7 > 3d/2 for the case described below in theorem[r|and 7 > d + 2 for theorem|2] 

We refer to the rys as random fields, although in general they may also be associated with some 
other parameters, e.g. random bond strengths as in the models studied in |[T1|3]. 

An example of a system of this type (with Kx = (Jx ) is the ferromagnetic transverse-field Ising 
model with a random longitudinal field ||T2| (henceforth QRFIM), with 

K = - E - E + + -i'^] (1-2) 

where a^^ {i = 1,2,3) are single-site Pauli matrices, and Jx-y > 0. When A = this is the 
classical Ising system with a random field. 

We will examine phase transitions where the order parameter is the volume average of the ex- 
pectation value of Kx with respect to an equilibrium state. These states are the infinite volume limits 
of sequences of Gibbs states in finite boxes K dl/^ with some boundary conditions, for which the 
expectation values are denoted here by {—) 

As is well known, this order parameter is related to the directional derivatives (±) of the (quenched) 
free energy density given by 

HT,Lh,e) = -T lim^log^A, (1.3) 

A^Z<i |A| 

where Za is the partition function of the system on a finite domain A of size |A|. Under fairly 
modest conditions ifTOl [131 , it is known that this limit exists, is independent of the sequence of 
boundary conditions and (almost certainly) of -q. Furthermore, it is a concave function of h, which 
means it has directional derivatives; and these are related to the order parameter by 

< hm ^ V (1.4) 



dh- - Ay-z" lAI ^ ' "^'^ - dh 
' ' xeA 

for any sequence of sufficiently regular (e.g. rectangular) domains A converging to Z**. 

More can be said, in particular in the classical case. There, convexity arguments imply that if 
the mean free energy density J-{T, h) is differentiable in h, i.e. 

dh^ ~ dh+ ^Im' ' ^ 
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then the block averages of over large rectangular volumes B do not asymptotically deviate from 

dh ■ 



the thermodynamic value — ^ . More explicitly, there is a positive function g{5) such that for every 



(5 > 0: 

>5\\ < Ae-g(*)l^l (1.6) 



lim Av ProbA 



1 v-^ dT 
2^^ 



iBi ^ " ah 

' ' xGB 



for some A < oo, where Av is an average over the random fields, and ProbA (— ) is the probability 
with respect to the Gibbs equilibrium measures for any sequence of A's which converges to Z"^. 
This general principle, which follows from the convexity of the free energy, can be used to show 
that the differentiability of the 'quenched' free energy density, i.e., ( |1.5| l which is the subject of our 
discussion, implies the vanishing of the 'Short Long Range Order parameter' as well as of the 'Long 
Long Range Order Parameter' which appears in ( |L4[ ). 

Combining the results known for some time for classical systems (where Hq and are functions 
of a local configuration) [.lOi . with those recently proven also for quantum spin systems (where they 
are self adjoint operators) ifTTl . we obtain: 

Theorem 1. In dimensions d < 2, any system of the form of ( |Ll| l has = for all h, and 
T > 0, provided e ^ 0. 



In view of equation ( |1.4| l, or the above comment, this means that no first-order (i.e. discontinu- 
ous) transition with as the order parameter is possible in the presence of randomness in the local 
value of h: any such transition present in the system described by Ha is "rounded". 

For systems where this spontaneous ordering would break a continuous symmetry, the results 
can be extended to d < 4. Here we must have observables which transform as vectors, so we 
replace with 

He = IC- ^{h + erj^) ■ ffx = "Hq - e^rj^ ■ ^x, (1-7) 

where h, fjx, and a-x are quantities which transform as vectors in a space of dimension N > 2 and 
K, is invariant under such rotations. 

Theorem 2. For the systems described above, with the random fields ffx having a rotation-invariant 
distribution, the free energy is continuously differentiable in h at h ^ Q whenever e ^ Q, d < A. 

As before, the proof for classical systems is found in |10| and the proof for quantum systems 
in mil. 



2 Higher Order Odd Correlations 

It can be proven rigorously for the models described by \\.2\ with purely ferromagnetic interactions 
(Jx-y > 0) and A = that in the pure (e — 0) system, not only the magnetization (^'^^^^ but all 

odd correlations, e.g. ^y^z, ai-e discontinuous at /i = for T <Tc IHl. The 

same is presumably true quite generally, particularly for the QRFIM with A 7^ 0; a discontinuity 
in the magnetization brings with it also a discontinuity in all odd correlation functions. We may 
now ask whether the randomness in the external field, corresponding to e 7^ 0, which removes 
the discontinuity in {ax) also removes the discontinuity in all the odd correlation functions when 
averaged with respect to the randomness. 
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Figure 1: Phase diagrams of (a) the simple Ising model and (b) Ising lattice gas in two dimensions; 
thick lines are first order phase transitions. 

For the cases covered by jl.l) when Jx-y > 0, this question can be answered in the affirma- 
tive, including for A 7^ 0. In this case the system satisfies the Fortuin-Kasteleyn-Ginibre (FKG) 
inequalities [15 1. Using these, it is possible to show that whenever all Gibbs states agree on the 
magnetization per site (a^) (i.e. when the spontaneous magnetization M is zero) they also agree on 
all other local observables (the argument for random systems is spelled out in 1 16 1). We do not have 
any such result for the case when Jx-y in ( |1.2| i contains also some negative terms. We expect of 
course that this will be the case more generally, but we lack a proof. 

For an instructive example let us consider the phase diagram of the Ising model ( |1.2| i with A = 
and d = 2 with nearest-neighbor interactions = J > (e e {(0, 1), (1, 0)}), e = (see 
Figure[T^). The free energy is analytic everywhere except on the bold part of the T axis. Along that 
line the free energy is continuous but its derivatives are discontinuous in all directions transversal to 
the T axis. That is, the free energy surface has a seam, i.e. a line of 'kink singularities'. 

Those singularities will be rounded through the addition of arbitrarily small randomness to the 
field h. However, the kink singularity would not be eliminated by small randomness in the bond 
strengths. That is, with J replaced by J + ebx^e with the bx.g i.i.d random variables, there would 
still be spontaneous magnetization at low temperatures, though T^. may change, and so will possibly 
the nature of the singularity in T atTc- The randomness is tangential to the seam, corresponding to 
"orthogonal randomness" in the terminology of |2|. 

Let us now consider this system in terms of lattice gas variables. Setting £^x = 5(1 + '^x) G 
{0, 1}, the nonrandom Hamiltonian now takes the form 

Ho = -4:J ^ - /I ^ + constant. (2.1) 

<a:,y> x 

The chemical potential /i is related to the parameters of the spin model by fi = 2h — 8 J, and so 
the phase diagram is now of the form shown in Figure [TJ?. At the bold line the free energy has a 
discontinuous derivative in both variables, and randomness in either J or /i will eliminate this first 
order transition, rounding the seam in ^{T, /i, 0). To avoid destroying the first order transition we 
have to add randomness to J and /i in a correlated manner, so as to make it 'tangential' to the seam, 
i.e., in the form done in ( |1.2[ ) in the cr-variables. 

Similar reasoning can be applied quite generally. In particular if we consider the free energy 
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Figure 2: Sections of the phase diagram for the Blume-Capel model. 

^(T, J,h,e = 0) for the system described by ( |1.2[ ) with J containing interaction parameters mul- 
tiplying both even and odd functions of the spins then, for T < Tc, T will have "multidimensional 
seams" in directions transverse to all the odd terms - which is why they are all discontinuous in 
the absence of randomness. Adding then randomness to any odd term, e.g. to the one spin term as 
in ( |1.2| i, then rounds the discontinuity in all of them. If on the contrary we add noise only to the 
even terms, this is tangential to the "seams" and doesn't destroy the first order transition. 

This is also the case when adding noise in the QRFIM ( |1.2| l to the transverse component of 
the field, i.e. letting A^ui^' — > + ^'r\x)<yx'^ ■ As proven in 1171 this does not destroy the 

discontinuity in (a!^''^. 

The situation is similar in the case when the spin at site x, S^, assumes the values = —I, 0, 1. 
A particular version of such a system is the Blume-Capel model, whose nonrandom version is de- 
fined by 

Ho ^ -J - ^y^^ - 5^ - /iJI (2-2) 

<x,y> 

The low-temperature phase diagram of this system is known rigorously from Pirogov-Sinai the- 
ory ifTslfTOl . A section is shown in Figure [2] alongside the corresponding diagram at zero tempera- 
ture. 

In analogy with the Ising case, randomizing the magnetic field h, hy h i-^ h + erjx, will destroy 
all the first order lines in d < 2 (in d = 1 this would be relevant only at T = 0). Randomness in J 
or A on the other hand should not destroy the discontinuity along the h = line. This is in accord 
with the analysis in [ 1 ]. Similar results will hold for the quantum version of this model with Sx now 
interpreted as si'^\ and a transverse field term J2x ^^^^ added. 

3 Effects of bond randomness on the discontinuity of sponta- 
neous magnetization 

While the 'tangential disorder', in the sense discussed above, will not eliminate the discontinuity in 
the order parameter associated with the transversal field, e.g. the spontaneous magnetization at low 
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enough temperatures, it may eliminate the possible discontinuity of the spontaneous magnetization 
if one was present at Tc in the non-random model. A general statement to that effect was proposed 
inim. 

An relevant example is provided by the two dimensional Potts models with nearest neighbor 
ferromagnetic interactions. At Q > 4 these models exhibit a discontinuity in the spontaneous 
magnetization at Tc, which coincides with a jump in the energy density. While Theorem 1 assures 
that randomization of the bond strengths will eliminate the energy discontinuity, it does not address 
the persistence or not of the discontinuity in the spontaneous magnetization. In [ 1 1 it was argued 
that also the latter discontinuity will be eliminated, but this remains an interesting open challenge 
from the rigorous point of view ll22l . 

4 One dimension, inverse square interaction 

Another class of models exhibiting a jump in the spontaneous magnetization are the one-dimensional 
Q-state Potts models, at Q > 1, with ferromagnetic pair interactions Jx-y which asymptotically 
behave as /\x — y\'^ . It is known that such systems have a phase transition, and that these exhibit 
somewhat unusual features ll20l 1211 (and references therein). Among those is the discontinuity of 
the spontaneous "magnetization" M(T) at zero "magnetic" field, which jumps at from zero to a 
value satisfying APj+/Tc > 1. 

At the moment it is not clear to us whether for the Potts models the discontinuity in M(T) 
may coincide with a discontinuity in the energy. That will of course correspond to a jump at Tc in 
the derivative of !F with respect to T. If that occurs then the system will have T — Teat least Q + 1 
phases. In the absence of energy discontinuity one may expect that at the critical temperature the 
systems will have Q ordered phases in analogy with the results of 1 14| for Ising systems. 

Let us consider now the effects of the additions of a random symmetry - breaking "magnetic" 
field and of randomization, at suitably decaying rates, of the spin-spin interaction terms. As men- 
tioned in the Introduction, the proof of rounding phenomenon |10| applies also to the long range 
d = 1 models as long as the pair interaction decays faster than \x — y|^3/2+(5 ^^j^jg ensures that the 
interaction energy between an interval of size L and the rest of the system is less than L^/^). Thus, 
in the presence of random magnetic field the spontaneous magnetization M{T) will vanish for all T. 
Randomization of the pair interaction would remove the discontinuity in the energy, if there was one, 
but unlike for the case discussed above we do not expect bond randomness to remove the disconti- 
nuity in M (T) at the (possibly altered) critical temperature. The reason for the possible difference 
with the two dimensional models is the bound Af^ J+/T > 1 which relates the discontinuity in the 
d = 1 models to the purely long range part of the interaction. 
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